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Abstract 

The plasma produced by a high-power pulsed circularly polarized microwave field which 
is much weaker than the atomic fields is investigated in the non-relativistic limit. It is shown 
that the resulting electron distribution function is non-equilibrium and anisotropic. Further- 
more, the dispersion equations derived to study and analyze low-frequency ion oscillations. 
It is shown an instability may grow in the aforementioned system only when the ion density 
exceeds a critical value. 
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I. INTRODUCTION 



A gas discharge in a microwave (MW) field is a rather comphcated phenomenon that exhibits 
a variety of features. Numerous theoretical and experimental studies have been devoted to this 
phenomenon. Interest in such a discharge is related not only to its potential for technological 
applications, but also for general understanding of the discharge physics. ^'^ A large amount of ex- 
perimental data has been accumulated and many theoretical studies have been made to explain its 
various details. A number of experiments have been performed with radiation sources operating 
at wave lengths of 8.5, 4, 2, and 0.8 cm.^ The strong MW field interaction with a gas raises inter- 
esting questions about the ionization mechanism and the electron distribution function created by 
this field. Furthermore, the interaction between high-power microwave fields with inhomogeneous 
plasmas results in a variety of interesting phenomena such as frequency upshift.^ 

In such a strong microwave wave field the electron oscillatory energy ee is much higher than 
the ionization energy /ioniz S^^ atoms 

where cuq is the radiation frequency, m is the electron mass, and Eq is the electric field strength. 
When the field amplitude is comparable to the atomic field strength £'a ~ 5.1 x 10^ V cm~^ 
tunneling ionization becomes an important mechanism for direct ionization of gas atoms. This 
effect has been thoroughly studied in previous papers. While condition (1) is satisfied by a 
large margin, relativistic effects are unimportant since the kinetic energy of electron oscillation 
in an electromagnetic field is comparable to the electron rest mass energy, for present day pulsed 
sources. In this case, gas atoms are primarily ionized by the impact of oscillating electrons and 
ionization is governed by the development of an electron avalanche. Recently it has been shown 
that discharge plasmas produced by super strong pulsed MW fields are subject to the Weibel 
instability associated with the anisotropic nature of the electron distribution function.^ 
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In the present paper, we consider the interaction of circularly polarized MW pulsed fields with 
frequencies a;o ~ 2 x 10^° — 2 x 10^^ s~^ with a neutral gas in the non-relativistic regime; in 
particular we study low frequency ion oscillations and their instabilities. We will consider a pulsed 
radiation source which is capable of generating radiation with an intensity of about 10^ W cm~^, 
whose electric field Eq < 10^ V cm~^ is much weaker than the atomic field E'a- In this case, we 
study the electron distribution function (EDF) and the stability of the discharge plasma in the 
aforementioned frequency range at non-relativistic electron oscillation energy ee . 

This work is organized in five sections. In Sec. II we review the EDF generated by the 
interaction of linearly polarized MW fields with a neutral gas. In Sec. Ill we consider the effect 
of circularly MW fields on the plasma produced and find the general dispersion equation for low 
frequency oscillations. In Sec. IV we solve the dispersion equations for low frequencies to find an 
instability of the system and its growth rate. Finally, a summary and conclusions are presented 
in Sec. V. 

II. DISTRIBUTION FUNCTION 

Under condition (1) the thermal velocity of the electrons in a discharge plasma can be neglected 
in comparison to the electron oscillation velocity in the MW radiation field. Since the collision 
frequency is much smaller than the MW field frequency, we can ignore the coUisional randomization 
of the forced electron oscillation as well. Furthermore, if the plasma density ne{t) produced by 
the field during gas breakdown is less than the critical density (i.e., uJq > uj^q = 47rne^/m) we can 
neglect the effect of the polarization field. Moreover, the plasma density is assumed to be much 
less than the neutral gas density no so that the latter can be considered constant. We also suppose 
that the field was adiabatically switched on in the infinite past. In addition, we can assume that 
the MW radiation electric field Eq is constant during a single field period. Therefore, by assuming 
that the electric field depends only on time, the kinetic equation for the plasma electrons produced 
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in the gas breakdown by a strong pulsed field can be written as follows^ 

^ + eEo • ^ = noa;ioniz5(p), (2) 

in the non-relativistic limit. Here 5(p) is the Dirac delta function, to'iQ^i2 the ionization probabil- 
ity of the gas atoms, and Uq is the neutral gas number density. In the case of MW gas breakdown, 
electron- impact ionization is governed by the ionization probability uj^q-^^^. From Eq. (2) the 
avalanche ionization parameter 7(£^o) can be determined We can neglect the right-hand side 
of Eq. (2) in the first approximation and calculate the EDF directly by solving the Vlasov equation 
under the following condition, 

^0 > 7(^o),^^ioniz- (3) 

Condition (3) depends strongly on the neutral gas density and is well satisfied at gas pressures 
of po — 10 — 100 Torr. In this approximation, to calculate the electron energy distribution function, 
we assume the field components with circular polarization to be 

Er^ — Eosincuot, Ey — Eq cos cuot, E^ — O, (4) 

where the electric field amplitude Eq describes a slowly varying (over the field period) MW pulsed 
envelope. 

Considering condition (3) and solving Eq. (2), by averaging the distribution function over the 
MW field period, we find the following expression for the EDF^ 

Mv±,v,))^ ^f'l (5) 

where, /o(v, t) = ne(t)/o(v) and ve = eEQ/imuo is the electron oscillatory velocity in an alter- 
nating electric field; Vz is the perpendicular component of electron velocity with respect to the 
electric field; v± is the component of velocity in the plane of the electric field; /o(v) satisfies the 
normalization condition / (iv/o(v) = 1. 
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Since a pulse of a strong MW field actually contains on the order of 10-100 wave periods the 
EDF is in principle unsteady. However, as one would naturally expect the form of the EDF is 
determined after several ionizatiug collisions. Thereafter the electron concentration exponentially 
increases in time while the shape of the EDF does not change. It is obvious that just on the 
initial moment of the gas breakdown the EDF is not in equilibrium. But after several periods 
an EDF is produced such that it usually characterizes the breakdown independent from initial 
conditions. Furthermore, the growth rate of the plasma density is expressed by ~ ^/ujq. Thus, 
under the gas low-pressure condition (3) the velocity distribution function of electrons is formed 
within the field period and is not practically changed furthermore. In addition, during instability 
development interested in the present paper the plasma density can be considered constant. The 
rate of collective effects grows with the increase of the plasma particles concentration, and this is 
why the growth rates of plasma instabilities, usually small compared to the plasma frequency, are 
proportional to the plasma frequency. 

Under conditions (1) and (3) the EDF is well described by Eq. (5) since, at the instant of 
originating, the electrons are evenly distributed in phase with the wave electric field. 

III. DISPERSION EQUATION 

The theory of the interaction of MW fields with a plasma has been widely developed and covers 
an increasing range of plasma phenomena.^^ In Sec. II the distribution function of the charged 
particles for a homogeneous plasma produced by an external circularly polarized MW electric field 
defined by Eq. (4) was obtained. The stability of such a plasma can be analyzed by considering 
small perturbations of this distribution. To begin with, we analyze MW fields with frequencies 
higher than all characteristic plasma frequencies 

<^0 > t^pa, l^a , (6) 
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where ut-pa, are the plasma frequency, and coUision frequency for a species (ions or electrons), 
respectively. To first-order approximation, the plasma in a MW field may be regarded as isotropic. 
Consequently, the field Eo(t) obeys the dispersion equation for transverse wave 



fco = — 
c 



;2 

'•pe 

ul C 



1 - ^ ^ -• (7) 



where c is the speed of light in vacuum. Therefore, the MW field can be considered homogeneous 
for processes with a characteristic length of inhomogeneity (where is much smaller than 
A;^^, or k ^ ojq/c). Assuming this condition, we analyze quasi- longitudinal oscillations in the 
stability of the MW-field-produced plasma. In reality, the oscillation frequencies are of the order 
of characteristic plasma frequencies, i.e., uj ~ a;pe. Taking into account relation (7), we find the 
validity condition of the quasi-longitudinal approximation to be a; <S a;o — k^c <S kc. 

The MW electric field with circular polarization primarily interacts with the neutral gas and 
produces a plasma whose EDF and properties were summarized in Sec. II. Since we will study 
instabilities whose growth rates are higher than the ionization rate we can assume a constant 
electron density. It is clear that, for relatively weak fields in real plasmas, electrons are accelerated 
to velocities significantly exceeding their thermal velocity during the period of oscillatory electric 
field and the mean free time. The result is a state in which electrons are moving with respect 
to motionless ions. Such a charged particle beam is unstable and its growth increment is rather 
large, of the order of the plasma ion frequency or even larger. In this section we derive the general 
dispersion relation for the MW produced plasma in the limit of low-frequency ion oscillation 
perturbations. Since the effect of a high frequency electric field on the plasma ions can be neglected 
compared to its effect on the electrons, the ion velocity distribution function can be regarded as 
isotropic (usually Maxwellian) and the EDF is given by Eq. (5). 

Linearizing the Vlasov equations for electrons and ions in a coUisionless plasma we obtain the 
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following equations 
dSfe 



, , xf <^^0r ■ +^ dSfe eB a/oe(p-Po) „ 

_ + • V d/e H (ismcuo^ + j cosa;o^) • -7{ \ 7{ = 

at m av m ov 



(8) 



(9) 



dt M dv 

for small deviations from the zero-order distributions. Here, i and j are unit vectors along x 

and y direction, respectively; /oe(p ~ Po) and /oi(p) are the zero-order distribution functions of 

the electrons and the ions, respectively; e , m, ej, and M are the electron and ion charge and 

mass, respectively and po = e J^'Eodt — — eEo cos ujot/ujo and i — \J—\. The small perturbations 

5/e and bj\ depend on the coordinates in the form of exp(^k.x), because of the assumption of 

homogeneity; E = — V0 is the electric field due to the perturbations satisfying the Poisson 

equation 

k^4> = 47re J dp Sfe + ^Tre^ j dp Sf^. (10) 
In order to solve Eqs. (8)-(10) it is convenient to introduce a new function 

V'e( p,t) ^ expl-^^^^ {kx sin uJot + ky cos uJot)\ /e(p + — [icosa;ot - jsina;oi]), (11) 
( muQ J ujo 

where k^ and ky are the component of the wave vector k. The system of equations then takes the 

form 



X 



dt 



dv 



e Jdp^e{p)+ei IdpSfi{p) exp 



mujr 



(kx sin ujot + ky cos uot) 



0, 



X 



dt 9v 

Cj / dp 5f[ + e J dp i/je exp \ ~ _ 2 sinujot + ky cosojot) 



'121 



mw, 







0. 



Applying the expansion 

exjp{±i{kxrE sincuot + kyrECOScuot)} 



J2p=-oo Em=-oo exp {±z{p + m)uJot} CXp|±Z — I JpikxTs) JmikyTE), 



(13) 



and introducing the decomposition 



|-oo +CO 



= ex.p{—iu!t) ^ ^ exp + n)u!ot} 

;=— oo n=—oo 



(14) 



we find from Eq. (12) the foUowing equations 



X 



-I [a; + (Z + n)u^\ + • v Ve,„ - «k • ^ ^ 

^ ^ Jl-p{kxrE) Jn-m{kyrE) 

df ■ Attp- 



0, 



X 



q / c^p(54^ + e /dp E+r-oo E+^-oo V'e.„ exp 



im — n)7i 



Jp-l{kxfE) Jm-nikyTE) 



= 0, 



(15) 

where r^; = eEo/ma;Q is the amphtude of the electron oscillation in the MW electric field and J„ 
is the Bessel function of order n. 
Now using the notation 



^&in jdp ipet„ , Ui^^ = ei j dp Sf^^^ , 



we can write the formal solution of Eq. (15) in the form of a coupled set of equation, as follows: 



ue,^ + Ejf-oo E^-oo exp 



'^ei^ = -^^e (a; + (/ + n)a;o, k) x 
(m — n)7r 



Jp-l{kxfE) Jjn-n{kyTE) 



(16) 



— —5ei (a; + (Z + n)uQ, k) x 



■"ii^ + Epz^^oo Em=-oo ''^ep^ exp "I — ^ ^ 



where fee (i^, k)andfei(c<j, k) are the partial contributions of the electron and ion susceptibilities 
to the longitudinal dielectric permittivity* 



£(a;,k) = l + fei + fee- 
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The existence of solutions to this infinite set of coupled equations gives the dispersion equation 
for small- amphtude longitudinal oscillations in the external MW electric field. Since it is a deter- 
minant of infinite order, it is impossible to analyze in general. However, for the most interesting 
limiting cases this determinant can be simplified and the oscillation spectrum can be found. For 
example, all the values Se^^lcu + ncuo, k) with n ^ 0, are small compared to unity in the very 
high frequency limit, when condition (6) is satisfied. Assuming this limit, only Ueoo and U[^^ are 
nonzero in Eq. (16); thus, we obtain 



[1 + dee{uj, k)] Ue^o + See{uj, k)Jo{k^rE) MkyVE) u^^^ = 0, 

JoikxTE) MkyVE) Se^{uj,k)ueoo + [1 + fei(a;,k)] u^^^ = 0, 
which gives the following dispersion equation 



(17) 



1 + deeiuj, k) + Seiicu, k) + dse{uJ, k) dsi{uj, k) [l - J^{k^rE) ^^( V^)] = 0- (18) 

By deriving Eq. (18), we fully neglected particle collisions. However, this equation is also valid 
for coUisional plasmas if the oscillation velocity component of the particles is small compared to 
the thermal velocity. The external MW field has no effect on the collision process, i.e., on the 
scattering cross sections. 



IV. INSTABILITY OF THE MW PRODUCED PLASMA. 

The instability analyzed in this section occurs in the range cuq ^ kvE ^ (^pe- It only occurs for a 
positive (non-equilibrium) slope in the EDF (5) and never for an equilibrium (Maxwellian) electron 
energy distribution in this frequency range. Taking a Maxwell distribution function for cold ions 
and the distribution function (5) for electrons, the electron and ion susceptibilities SeQ[{u!, k) 
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are given by*'^^ 



. . ,^ "^pe f, kdfo{v,t)/dv uj^ek 



?>.A;2 s ool^s (19) 



Lip V 



where 



1 

1/2 



for 
for 



i?e(w) 7^ 
it:e(u;) = 



when 



E 



< 1 



4k^,v^ 



> when 



> 1 



for Im{ijj) — 

1 for 7m(a;) 0, i?e(a;) 
1/2 for i?e(w) = 

where Re{uj) and Im{ijj) are the real and imaginary part of a;, respectively. Substituting expres- 
sions (19) and (20) into the dispersion equation (18), we obtain the following dispersion relations 



(20) 



1 + ^ 



4:k]_v% > co- 



pe 



\2k±VE\ k 



± 



k 



( \2k^VE\' k ^pi 



1-1, 

2 \^ a;2g k 



kj_ 



corresponding to unstable low frequency ion oscillations. Here, the following approximation 



r 



{21) 



(22) 



was used. Equation (21) show that in the case of a circularly polarized MW field the low frequency 
oscillation takes place only for perturbations satisfying the following condition: 



UJ. 



\2k^VE\^ ~ A;2 



pi 



(23) 



This instability resulting from the positive slope of the distribution function (5) can be interpreted 
as the stimulated Cherenkov excitation of the low frequency ion oscillations by the radiation field 
of the oscillating electrons. Perturbations that do not satisfy the condition (23) are not unstable. 
Of course, according to approximations and assumption used, the instability satisfying condition 
(23) can only occur when the growth rate exceeds the ionization rate, i.e., the growth rate should 
exceed the avalanche ionization constant 7(-Eo)- 

V. CONCLUSION 

In analyzing the interaction of a high frequency circularly polarized MW field with a neutral 
gas, we obtained the general dispersion relation for low frequency waves in the MW produced 
plasma. Prom Eq. (23) for the instability condition it is necessary that the the ion density exceeds 
the critical value Ui > Mk\EQ/T{w?k'^ujQ. Under such condition the electrons transfer electric 
field energy to the ions and stimulated Cherenkov excitation of low frequency ion oscillations 
takes place. This effect requires a positive slope in the distribution function and, in this case, the 
instability can propagate when the grow rate exceeds the avalanche ionization rate '^{Eq). Other 
work underway, is concerned with the case of linearly polarizd MW fields. 
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